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We demonstrate the simultaneous quantum state reconstruction of the spectral modes of the 
light field emitted by a continuous wave degenerate optical parametric amplifier. The scheme is 
t^j- \ based on broadband measurement of the quantum fluctuations of the electric field quadratures and 
subsequent Fourier decomposition into spectral intervals. Applying the standard reconstruction 
algorithms to each bandwidth-limited quantum trajectory, a "spectrum" of density matrices and 



Wigner functions is obtained. The recorded states show a smooth transition from the squeezed 
vacuum to a vacuum state. In the time domain we evaluated the first order correlation function 
of the squeezed output field, showing good agreement with theory. 



The experimental techniques of quantum state reconstruction, first applied five years ago, have 
c5 ' opened a new field of research, wherein simple quantum mechanical systems can be characterized 
completely by density matrices and Wigner functions Our system consists of electromagnetic 
field modes at optical frequencies. We have previously generated the whole family of squeezed 
states of light using an optical parametric amplifier (OPA) and reconstructed these states using the 
method of optical homodyne tomography ||. The reconstructions presented therein were limited 
to essentially one particular pair of modes at frequencies where Q is a radio frequency and u 

the optical frequency. The spectral bandwidth of this mode pair was AQ/2ir = 100 kHz. Since an 
OPA pumped below threshold emits a frequency spectrum, with a bandwidth determined by the 
cavity linewidth in the order of several MHz, the output of the OPA is described more precisely by a 
whole spectrum of quantum states [0] . General schemes for multimode reconstruction become quite 
complicated already at the two-mode level M. So far only one experiment demonstrating photon- 



number correlations in the time domain by measuring the photon statistics via dual-pulse phase- 
averaged homodyne detection has been carried out 0. In this work we use a simple measurement 
scheme to record simultaneously the complete statistical information about the quantum states of 
all mode-pairs emitted by the OPA, disregarding correlations between the different mode-pairs. 
This enables us to reconstruct the Wigner functions and density matrices corresponding to a large 
set of mode pairs. 

The experimental setup is shown in Fig.|l| Center piece of the experiment is a monolithic 
standing-wave lithium-niobate OPA, pumped by a frequency-doubled continuous- wave Nd:YAG 
laser (532 nm). The output of the OPA is analyzed by a homodyne detector whose difference 
current z_ is directly proportional to the OPA's output electric field. In order to obtain simulta- 
neously information about the quantum states of all modes emitted by the OPA, the homodyne 
detector current £_ was recorded with a bandwidth covering the whole frequency range up to 30 
MHz, exceeding the OPA's 17.5 MHz cavity linewidth. For this purpose fast photodiodes (Epitaxx 
ETX500T, specified bandwidth 140 MHz) with broadband amplifiers and a 30 MHz A/D-board 
(IMTEC) for data collection were employed. The z_ data (about 500 000 points with 12 bit res- 
olution) are taken while the local oscillator phase is swept by 2tt in approximately 8 ms. The 
recorded noise trace is split into 16 separate traces, each containing the information of a mode pair 
of small bandwidth 1.9 MHz, centered at a sideband frequency u±Q. This is done by dividing the 
Fourier transform of the recorded trace into 16 equal intervals, and taking the inverse transform 
of each interval separately. The same is done for a measured i_-trace of a vacuum state, in order 
to normalize each of the 16 traces. Due to electronic excess noise of the detection system at low 
frequencies the first of the 16 traces, containing the spectrum between and 1.9 MHz, is discarded. 
The quantum states of the remaining 15 modes can be reconstructed the same way as described in 
0, employing the inverse Radon transform ]IJ and the pattern functions of the harmonic oscillator 
0. Thus a whole "spectrum" of Wigner functions and density matrices is obtained from a single 
temporal record. 

The main difference between the squeezed vacuum states reconstructed at different Q is the 
amount of noise reduction (squeezing) and enhancement (antisqueezing). The spectra of quantum 
noise power of the squeezed quadrature and of the antisqueezed quadrature of an OPA on 
resonance are given by 



where d = ^ P/ P t h is the pump parameter with pump power P and a threshold power P t t of the 
OPA, T/27T = 17.5 MHz is the cavity linewidth (HWHM), ty is the spectral density of the vacuum 
state, £ = 0.88 is the escape efficiency of the resonator, and rj is the detection efficiency. The latter 
includes propagation losses after the OPA, homodyne efficiency and detector quantum efficiency In 
contrast to our previous work we did not employ a mode cleaning cavity, thus classical noise of the 
pump wave was not negligible and the modematching between local oscillator and the OPA output 
signal was ~ 95%. Furthermore the efficiency of the photo detectors had degraded slightly from 
the previously measured 95-97%, thus the overall detection efficiency including escape efficiency 
£ ■ rj was fitted to be 0.7. The pump power was approximately half the threshold power. Fig.|2] 
shows the spectra of the maximally squeezed and anti-squeezed quadratures measured directly 
with a spectrum analyzer in comparison with the ones obtained via multimode reconstruction. 

Three of the reconstructed states are plotted in Fig.[| They display strongly elliptical Wigner 
functions and oscillations in the photon number |7|], characteristics that vanish when the state 
aproaches a vacuum state at frequencies sufficiently away from the OPA cavity center frequency. 
We remark that we did not include error corrections in the reconstruction algorithm, as described 
in || . For a single mode treatment of previously measured data of us that takes into account the 
finite detection efficiency see |§ . Note also that the present measurement method does not allow to 
make any statements about possible correlations between the different mode-pairs. Regarding the 
light field generated by the OPA this is of minor importance, since according to theory mode-pairs 
at different frequencies should be uncorrelated ||. 

A direct measure of the total noise of the squeezed and anti-squeezed quadratures / \l/±(fi, P)dVL 
of the overall OPA output within the bandwidth 1.9-30 MHz is obtained by analyzing the full 
quantum noise instead of its particular Fourier components. (For this, the non-uniform spectral 
response of the A/D-board has to be eliminated via normalization in Fourier space.) The inset of 
Fig|| shows the total electric field variances as a function of the local oscillator phase. Well-known 
figures of this type usually depicted only the .E-field variances of a single mode of the light field [[ll]]. 
Here we obtain for the overall OPA output a total squeezing of -2.9 dB (= 0.47 linear scale) and 6.7 
dB (= 4.7) for the total antisqueezing. The corresponding photon statistics is shown as well. This 



multimode light field would be detected if one would employ photon counters instead of homodyne 
detection. Note that no photon number oscillations are observable in the total OPA output. This 
may seem surprising, since almost all the photon number distributions of the individual modes 
show oscillations and all the reconstruction transformations are linear. However, according to 
basic laws of probability theory the overall sampled distribution of the quantum noise is not given 
by the average but by the convolution of the distributions of the individual modes. Therefore the 
photon statistics of the total OPA output is not given by averaging the single mode statistics. 
What is averaged though in our case of uncorrelated Gaussian noise distributions is, as mentioned 
above, the variance at each particular phase angle. 

A further perspective is gained by the analysis of our data in the time domain. At a fixed phase 
9 of the local oscillator, the recorded time trace can be regarded as the quantum trajectory of the 
quadrature xg(t) = 2~ l l 2 (a) '(£) e id + a(t) e~ l9 ), where a(t) is the Heisenberg output field operator. 
Thus the recorded data contain all information needed to extract the first order correlation function. 
This function is easily calculated using the input-output formalism for optical cavities |§. 
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Here <>e means averaging over all phase angles. Fig.[| shows that the evaluation of the 
first order correlation function obtained from the experimental data is in good agreement with 
theory. In principle it is possible to obtain via broadband homodyne detection all higher order 



time correlations of the field quadratures {xe{Ti)xQ{T2)---Xe{r n )) . Except for one special case |TT 
their significance does not appear to have been studied. For our data the evaluation of these 
correlations could not be done reliably, due to difficulties with a proper normalization at r = 
and since the minimum recordable timestep of 1.7 ns was too large in comparison to the expected 
rapid exponential decay times e _n ( 1_d ) rr . 

Note also that since a homodyne measurement is always a two-mode measurement of the pair 
of modes at frequencies our experiment shows two-mode time correlations. For single-mode 



time correlation measurements a heterodyne system (detection efficiency rj <0.5) can be employed. 
This is possible, since detection efficiency does not play the same crucial role in time correlation 
measurements that it does for quantum state reconstructions, in fact ideally g^ n ) is independent of 

T). 

In conclusion, we have demonstrated the first simultaneous reconstruction of a whole spectrum 
of quantum states of the optical light field. Some of the reconstructed states emitted from the OPA 
display photon number oscillations and strongly elliptical Wigner functions. These characteristics 
vanish for the mode-pairs at frequencies sufficiently away from the OPA cavity center frequency. 
Furthermore the first order time correlation function was evaluated in good agreement with theory. 

Our measurement scheme may be also useful for a variety of other quantum optical systems 
with more complex frequency or time dependencies. Examples are recent squeezing experiments 



with solitons in a fibre ||12|| , pump-noise-suppressed laser diodes ||T3[, or exciton-polariton systems 
in semi-conductors []14,15|. 
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FIG. 1. Experimental scheme for measuring the quantum states of the OPA. 
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FIG. 2. Spectrum of squeezing and anti-squeezing: Spectrum analyzer traces (upper and lower grey 
lines) in comparison with a theoretical fit (smooth black lines) using eq.l. Each of the 15 pairs of 
black squares corresponds to one reconstructed quantum state at a specific frequency Q and represents 
the maximum/minimum width of the sampled marginal distributions. The reduction in squeezing at 
frequencies below 4.5 MHz is due to classical laser noise at low frequencies. 
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FIG. 3. Three examples of the simultaneously recorded squeezed vacuum states emitted by the OPA: 
Upper row: contours of the Wigner functions, lower row: corresponding reconstructed photon number 
distributions (dots) with theoretical expectations (histograms). For increasing offset frequency from 
the OPA cavity center frequency the states aproach the vacuum state and the characteristics of squeezing 
such as photon number oscillations and ellipticity of the Wigner function vanish. 
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FIG. 4. Electric field variances in dependence of the local oscillator phase and photon statistics of the 
sum of all OPA modes between 2 and 30 MHz. The average photon number (n) = 0.8 gives the average 
photon flux per Hz bandwidth. This implies a total photon flux of 2.2 • 10 8 photons/s « 45 pW within 
the detection bandwidth of 28 MHz. 
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FIG. 5. First-order time correlation function of the recorded noise trace (dots) in comparison with 
theory (line). A quasi-exponential decay similar to the one of a thermal state is observed. 
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